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I INTRODUCTION 


In the early 1970s, stealth technology was revealed to the public. Since then, the 
term stealth has generally been interpreted to mean “invisible to radar.” In reality, however, 
stealthy targets are not completely invisible. Stealth technology merely reduces the radar 
cross section (RCS) of the target. The RCS of a target is denoted by o and may be defined 
as: 

_ Power reflected to receiver per unit solid angle 


Incident power density / 41 


A great deal of effort is spent to reduce RCS using four basic approaches: 


e Target Shaping 
e Surface Material Selection and Coatings 
e Active Cancellation 


e Passive Cancellation 


Each of these methods involves trade-offs. For instance, the optimum target shape to 
achieve reduced RCS may result in aerodynamic problems: Additionally, RCS reduction 
methods (except for shaping ) are usually very narrowband and are effective over limited 
Spatial regions. These factors dictate that RCS reduction methods must be chosen based on 
the platform’s mission and potential threats. The main goal of RCS reduction is to achieve 
an RCS value below an acceptable threshold over a specified range of frequencies and 
angles. 

The relationship between the electrical characteristics of the material employed and 
the target’s resulting scattered field is a very important concept in RCS reduction (RCSR). 
The goal of materials selection is to choose a material that yields the smallest reflection of 


radar waves. The use of radar absorbing materials (RAM) is crucial to RCSR, but because 


of the complexity of the problem, RCS reduction using RAM has historically been based on 
intuition and experience. 

This thesis examines a method of synthesizing a target’s surface material’s electrical 
properties to obtain a specified RCS. Synthesis is a special case of the more general 
problem of inverse scattering, whereby information on an unknown target is derived from 
its scattered field. In RCS synthesis, the target shape is usually known; determining the 
surface material is the primary goal. The synthesis problem can be solved in either the time 
or frequency domains. This thesis addresses the problem in the frequency domain. 

The solution of integral equations provides a basis for obtaining surface impedance 
information from a target’s RCS. In general, both electric and magnetic currents can exist 
on the surface. Integral equations can be solved by using the method of moments (MM), 
which is a frequency domain technique. The MM solution employed here requires that the 
entire structure to be modeled be decomposed into a number of surface patches 
(subdomains), which must be small compared to the wavelength. There are two approaches 
to deriving the integral equations for surface impedance synthesis. These are the 
approximate method and the rigorous method. In the approximate method, the unknown 
surface impedance is expanded in a series of basis functions with unknown coefficients, and 
an estimate of the surface current is provided. The rigorous method requires that the 
coefficients of series expansions for both the impedance and current be determined 
simultaneously. 

Previous works on RCS synthesis have discussed the problem in terms of the 
approximate method [1 and 2]. The rigorous method 1s the focus of this thesis. In general, 
for the rigorous method, a nonlinear set of simultaneous equations must be solved. The 
present work will be limited to a resistive surface and the target geometry to a flat 
rectangular plate. Additional boundary conditions will be utilized to solve the synthesis 
equations for both sets of unknown coefficients. In this limited case the resulting system of 


equations 1s linear. 


Chapter II presents theoretical background required in the formulation and 
evaluation of the problem. Reflection and transmission coefficients, boundary conditions, 
and surface impedance concepts are explained. In Chapter III, the RCS of a rectangular 
resistive sheet is calculated with MM and the RCS synthesis procedure for the ngorous 
method is developed. Appropriate equations are derived and expressed in matrix form for 
the bistatic case. In Chapter IV, computer implementations of these synthesis equations are 
discussed. Chapter V summarizes the results of the synthesis and presents conclusions and 


recommendations for future research. 


Il. THEORETICAL BACKGROUND 


A. SURFACE RESISTIVITY 


The goal of radar cross section reduction is to eliminate, or at least attenuate 
significant scattering sources on the target. Thin layers of lossy materials are of special 
interest in RCSR design, primarily in the treatment of surface waves. Discontinuity 
boundary conditions are used to mathematically represent thin films. They are the resistive 
sheet and conductive sheet boundary conditions. The boundary conditions define the 
relationship between the current and the tangential field components on the surface. A 
resistive sheet, which is the primary focus of this thesis, is an infinitely thin imperfect 


electric conductor. Its conductivity is finite, which means o, < « , and it does not support a 


magnetic current (J ms = 0). For resistivity, the unit ohms per square (Q2/ sq ) is generally 
used. This unit is derived from the basic definition of resistance for the special case of 
length equal to width (2 = w ) 

ei 


R= = . 
> o.tw ot ) 





where o, is the conductivity of the material per meter and ¢ is the thickness of the 
material. 


There are two limiting values for surface resistance. The first is R,=0, which represents a 


perfect electric conductor. The second is R,=oo, which represents a surface matched to free 


space. 
The boundary conditions may be expressed in general vector form as [1] 
nx E(+)-nvx E(-)=0, (2) 
nx A(+)-ax A(-)=J, , (3) 
with: 


fix {ax E(+)}=-RJ, , (4) 
where 7 is the unit vector normal to the sheet, and is directed to the upper side of the sheet, 


and J.is the total electric current supported. Plus and minus signs represent upper and 
lower faces of the sheet, respectively. 

The electromagnetic dual of the resistive sheet is the magnetically conductive sheet, 
which supports only magnetic currents. A combined sheet consisting of resistive and 
conductive sheets is equivalent to an impedance surface, which is a surface where 
permittivity and permeability both differ from the surrounding medium [1]. The surface 
impedance approximation, also known as the Leontovich boundary condition, is: 


—_tan_ 


s = A’ (5) 


where E,, and H,,,are the total tangential fields on the surface, and7,is the surface 
impedance in ohms [3]. Solution of a scattering problem for the resistive sheet gives the 


solution for the conductive case by way of the duality principle. 
B. REFLECTION COEFFICIENT 


When a traveling wave reaches an interface between two different regions, it will be 
reflected or transmitted with the relative magnitudes of the two components determined by 
the impedances of the two regions. The reflection coefficient for a resistive thin film can be 
derived from the boundary conditions in equations (2), (3), and (4). The geometry is shown 
in Figure 1. Bistatic scattering for parallel and perpendicular polarizations are investigated 
separately. Specular reflection can be assumed _ because the sheet is infinitely thin. 

a: Parallel Polarization 


For parallel polarization (TM case), the incident wave is 
E, =0Eje FFT (6) 


where k is the incident wave vector. After some mathematical manipulations, the reflection 


coefficient is found to be: [1] 


aeleCOsu: 


" 2R,+7, cos0, 7) 


Diy 
The surrounding medium is assumed to be free space with impedance n,. 
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Figure 1. Reflection and Transmission from a Resistive Sheet 
oe Perpendicular Polarization 


For perpendicular polarization (TE case) the incident wave 1s: 
E, = § Elbe 6) 
and the reflection coefficient for the resistive sheet is [1] 


abs 


Sea oe 2 
DRECOSG 21), ©) 


TE 


CG. RADIATION INTEGRALS 

The radiation integrals are integral solutions to Maxwell’s equations, and are also 
called Stratton-Chu integrals. They can be derived directly by taking the curl of 
Maxwell’s first two equations, using the vector of Green’s theorem, and then integrating 
[3]. E- and H-fields are determined from current distributions by using the radiation 
integrals. On a resistive sheet the magnetic currents are zero. The field equations are 
further simplified when only the far zone fields are of interest. 

In the far zone (kr >> 1) and the E- and H-field vectors are orthogonal to each 
other and the direction of propagation resulting in TEM mode fields. The equations used 
in this section are derived by Balanis [4]. A rectangular coordinate system is used in 
formulating the solution for the rectangular sheet problem. The geometry is shown in 


Figure 2. 








xX 


Figure 2. Coordinate System for Scattering from a Planar Rectangular Plate 


The far-field spherical components of the E-field are: 








£.=0 
-1..—Jjkr 
jke~/ 
Eo =- N 
8 Atr No 6 > (10) 
and 
-7._—]kr 
jke™/ 
E,=- N 
d ee NolV 6 (11) 
where 
Nie II lJ, cos8cos$+J, cos@sin ] eo ideiay, « (12) 
N, = II - J. sind+J, cos¢| Peachey. (13) 
and 
r'cosy = x’sinOcoso + y’sinOsing = x'ut+ y'v=g. (14) 


The direction cosines are 
u=sinOcoso , 
v=sinOsing , 
and 
w=cos®. 

In order to find scattered fields, the current must be known. The method of 
moments will be used to find the current. In this approach, the unknown quantity is 
represented by a series with unknown expansion coefficients. The boundary conditions are 
imposed to obtain an equation for the current. After the current is found the scattered fields 
can be obtained from equations (10) and (11). This is referred to as the direct problem, 1.e., 
finding the fields from the currents. Determining the currents from the far field is the 


inverse or synthesis problem. 


D. RCS SYNTHESIS APPROACH 


RCS synthesis is a part of a more general problem of inverse scattering. In an 
inverse scattering problem, the scattering field is known and the objective is to extract 
current and resistivity information about the target from the scattered field. There are two 
types of inverse scattering problems: general and restricted. In the general case nothing is 
known about the target whereas in the restricted case some characteristics of the target are 
known. The radar cross section synthesis problem falls into the latter case. The target shape 
is known, but the surface materials and current are to be determined. Synthesis problems 
can be solved in the time or frequency domains. This thesis deals with synthesis problems 
in the frequency domain. 


At an observation point P, the scattered electric field will be [3] 











a | zs Z a View 
E(P) = —jkn, (I]1.0+% i, Sa Gan iP voli (15) 
where 
7 JAR 
G= Se 
ae and R 7 : | 


Primed quantities are associated with the source point; unprimed quantities refer to the 
observation point. In this equation, £,(P) 1s known, but the current and surface 
impedance are both unknown. This equation can be solved approximately or rigorously. 
In the first case, the current can be approximated and the equation can be solved for 
impedance. Equation (15) is sufficient to completely describe the problem. In the second 
case, both current and impedance must be determined. This requires additional 


information provided by the following equation: [3] 


10 


=e PPP AC) 22. Har) x lv « [[7,Gas 


+ jk J [2,(@' x J,)Gds' += 2 {fv \Z,f J PGde' (16) 
feiss a ~IM(V"-F,)vGas! —V x [[z.(4 <I Gas 


Equations (15) and (16) can be solved simultaneously for the current and impedance. 
E. METHOD OF MOMENTS 


The method of moments is a frequency domain technique. It is used for solving 
complex integral equations by reducing them to a system of linear equations. The MM 
technique requires that the entire structure to be modeled be broken down into subdomains. 
The surface is subdivided to a number of patches with dimensions small compared to the 
wavelength. Moment methods employ a technique known as the “method of weighted 
residuals.” Actually, the terms method-of-moments and method-of-weighted-residuals are 
synonymous [5]. 

The MM technique is popular because it can be applied to arbitrary bodies. 
Furthermore, integral equations are reduced to a set of linear equations, which are 
conveniently solved by using matrix algebra. The size of the matrix is directly related to 
the electrical size of the body. Large bodies result in large matrices, which affects 
computation time, and limits the size of the targets that can be handled in practice. 

The first step in MM for the direct problem is to represent the current by a series 
with unknown expansion coefficients. For a rectangular resistive plate in the xy plane, one 


possible representation for the current is: 


Coy > Kile 1732) (17) 
i 
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The vectors J* and J» are the expansion or the basis functions, which are complex in 
general. Selection of the basis functions is very important. They should be 
mathematically convenient — in other words, easy to integrate and differentiate. They 
should also be consistent with the behavior of the current to provide fast convergence 
[3]. The expansion scheme for the currents can be seen in Figure 3. Two basis functions 
span each rectangular subsection: one x-directed and a second y-directed. These two 
orthogonal basis functions are sufficient to uniquely define all possible current vectors on 


the patch. 





Figure 3. Three-dimensional View of Pulse Basis functions for the Currents J, and J : 
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Il. FORMULATION AND SOLUTION 


A. RCS OF A RESISTIVE SHEET USING THE METHOD OF MOMENTS 


In this section, the scattered field and radar cross section of a resistive sheet is 
calculated for the bistatic case. This direct problem solution is used to generate and 
validate the synthesis procedure. A constant resistivity is assumed on each patch, and the 
method of moments is used to find the current distribution on the plate. Rectangular patches 
and nonoverlapping pulse basis functions are used for calculations. After the synthesis 
procedure, the resulting resistivity is compared with the original resistivity to verify the 
synthesis equations. 

The resistive sheet lies in the xy plane. The incident wave is assumed to be a planar 


wave, and it contains both parallel and perpendicular components [3], 
mes A : Ak es net nee es 
a (6 Ei +62, ]e jkr (2B + SE’ Je ay (18) 
where, 
—jk-7= jk|xsin®, cos, + ysin@,; sing, +zcos8, | (19) 
and, 


ES = E4cos®. cosh; — E} sind, 
ie i i b i 
20 
E! = Ei cos sing. + E; cos = 
Need i i ~$ i 


The resistive sheet integral equation is solved to find the current distribution on the surface 


[3], 


E,(F)) = RJ, +4 jo | J,Gds'+ Fv Ve [[I,Gas' (21) 
tan Ss 0 S tan 


IS 


The unknown current is expanded into a series using equation (17), which can be written in 


the following manner: 


J, = Ja yy = me oe] + Da ae (22) 


The quantities Ax and Ay are defined in Figure 4, and p,, ee y) is the pulse function 


] 
ey . ‘o 


The method of moments is initiated by inserting equation (22) into (21). At this 


, if x' and y' on patch n. 
, else 


point a set of testing or weighting functions are defined. The weighting functions describe 
how the boundary conditions are enforced [6]. These functions are generally chosen to be 


the complex conjugate of the basis functions ( W = J : ), which is known as “Galerkin’s 
Method” [1]. Weighting functions for this case are [7], 


W =W,2+W,y (24) 
where, 

N 

“i. aie ceny) 

W - 

f a a (25) 
N 

SR) 

Wy, = 
yD. . (26) 


Both sides of the equation (21) are multiplied with these weighting functions and integrated. 
In terms of expansion and weighting functions, the integral equation (equation 21) takes the 


form 


[[, e Wds = [[RJs e Wds + jkno [fas fas 3, 7 -L(veMv'eJ,) (27) 


Equation (27) can be separated into two equations, one associated with the x 
components and the second with the y components. The derivatives of the pulse functions 


are approximated by 6 functions as shown in Figure 4. The x component of equation (27) 


becomes, 
S(x -x,, y) 
S(x, y- y,) 
My 
oa 
D, 
A 
Yn =Vn ZF 
— (x, y— yz) 
Xx 


~5(x—x*, y) 


Figure 4. Pulse Function Derivatives for a Patch [7]. 
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A similar equation can be obtained by testing with the y-directed weighting functions, 


i Eye JE on . i al P pale) _—s - 


oe (jas far Pn 2) om x ale {a(x.9 - yz) 


-(s.7- 5] foes) -afe-at ro“ err 
-a(2',y’- v2) oe, yy). 


Now define the following quantities, 


Gan = 


= © [[dedty [[ae'dy'G(x,x',y.y’) 
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Ox = alt fay [ay \C(xn xn -y.")-G(xnoxt yy’) 
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Aym Ayn 


-G(x3.%; Voy ‘)+ G(xs,.xt.».y'} 
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(30) 


(31) 


(32) 


(33) 


(34) 


(35) 


(36) 


as Peer Jer"{olsz.vars")- Olas, yan) 7 


-G(x,x; Ynoy')+ G(x.x7 ny"! 





Gra 
g2 a (38) 
0 ,limen 
Va = B4 aysine{ — )sinc{ AAP dCs mi) (39) 


Using equations (30)-(39), equation (28) and equation (29) can be cast in matrix form : 


Vv : g° lV’ G* +Q™ Q” I" 
Veale 3g lV o Q” GY +Q” | P (40) 


For the synthesis problem, the number of unknowns in this equation is 4 N where JN is the 


number of expansion terms. There are 2N current coefficients and 2N __ resistivity 
coefficients in general. If the material is isotropic then 7° = 7” and the number of unknown 


resistivity coefficients is just N . For the direct problem the number of unknowns is only 
Ne 
The integrals appearing in the above equations must be evaluated numerically. To 


reduce the computer run time, the following approximation can be applied [3], 


A 


2 
[F(x.x')de = A ee!) (41) 
A 


Xm = 


Z 
With this approximation, the impedance matrix elements in equations (31), (32), (36) and 
(37) are 
_~ ~JNo 


Onn = z 1G(xinsn Ym Yn) — (m8 YoneYn) 


(42) 
—G( i .% 2 Ym) + G( x3, YmsIn)} 
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eel - - - 
Onn = ; 1G (Xin tn Ym Ve |- G(XinXn Yaa 


(43) 
—G( x5 5% ns YmneYn)+ G(X %nsYms de} 
_ JN — _ 
Os — eG Xn, YinzYn)— G(XmsXn YneYn] 
(44) 
—G(XnoXn YineIn) + G(x x2 YnsYn)} 
al a - 
Q:. = ENO Gx ns Ye Yn )— OS mo%nsYins Ii) 
(45) 


-G(Xms%nsYarYn | + G(Xins%n Ymva) 
Note that in all of these equations, the Green’s functions in the integrands can 
become singular for the self term because both the source and the test points lie on the same 
segment. This problem can be solved by moving the test point off the patch axis [3]. For 

example, when integrating in x or x‘ 

he 1\2 ’ - i 
= (x—- x’) +(y-y’] uae (46) 
and when integrating in y or y 
2 

R= Gaa) a=, ) +a,” (47) 


In this case, an equivalent radius of the patch is used : 
a, = 0.225 Ax G0 225 hy (48) 


Equation (40) provides a means of solving the direct problem by plugging the 
resulting current coefficients into the far-field radiation integrals. The current coefficients 


can be obtained if the incident field and resistivity are known. The scattered field at an 


observation point (0 529 es) is computed by integrating the current. The current coefficients 


19 


that occur in the scattering equations are the same for all observation angles, because the 
same incident field determines these coefficients. After integrating, the scattered electric 


field has the following form: 


— , a N ° 
E. = ~JIMo ike [2 Ty Ax sind SA) sin 2) 





Amr z 2 
n=] 
Ree . ( kAxu,\ . ( kAw k 
+ ~) IJ Ay sine : - sine{ HAs a (49) 
n= 


= XE, + YE, 
where: 
§ = Xpls TVnVs 
= x, sin@, cosd, + y, sin®, sing, 
Finally, the complete expression for the 8- and - components of the scattered electric field 
IS: 


Ea = E, cos0, cosh, + Ey, cos®, sing, (50) 
Ey, =—E, sino, +E; cosd, (51) 
From the definition of the radar cross section fora three-dimensional target [1], 


ae 
Z z, 

o = lim| 4z2r 7 ape 
ro E; 


(52) 


and the corresponding radar cross section for the 9- and o- components of the scattered 


electric field for a unit plane wave are 





kn] » s]? 

S99 =| Bo" (53) 
kno | 5]? 

90 ag Le] 64) 
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The first subscript on o denotes the polarization of the scattered field, whereas the second 


one denotes the polarization of the incident wave. 
B. SYNTHESIS FORM OF THE SCATTERING EQUATIONS 


The resistivity of the rectangular plate can be determined rigorously by expanding 
both current and resistivity into series and solving for the expansion coefficients. For this 
purpose, the scattered electric field is found using the formulation in the previous section. 
Because the sinc functions in equations (34) and (39) are periodic in direction cosine space 
(DCS), it is convenient to work in DCS (Figure 5). There are two approaches sampling the 
RCS in DCS. The first is to choose the observation points entirely inside of the unit circle. 
The second is to include observation points outside of the unit circle. The first approach 
makes physical sense because all observation directions are points in real space. Directions 
outside of the unit circle correspond to complex angles. However, the first method 1s 
restrictive in the number of observation points, because they must be equal to the number of 
subdomains on the plate to obtain a square matrix that can be solved by inversion. The 
computer program which calculates and plots RCS in DCS 1s presented in the Appendix. 

In synthesizing the plate resistivity the far-field pattern is given. It is known that a 
resistive sheet does not support magnetic currents, and the far-field is related to the electric 
currents on the sheet by 

E,(P) = —jkmo || J,Gas" (55) 
s 
P denotes an observation point in the direction (r,0,6) or (oa, y, z) . On the surface of the 


sheet, equation (21) still holds, and in this equation resistivity and current are not 

independent of each other. These two quantities are related at every point on the surface of 

the sheet by the tangential component of the electric field as indicated in equation (4). 
Equations (21) and (55) can be solved simultaneously for the current and resistivity. 


For the current expansion equation (22) is used, and for the resistivity, 


a 


N 
Ry = > RaPnlX'ny’). (56) 


n 


Unit Sphere 
( r=1 ) 





Figure 5. Direction Cosine Space (DCS) . 


The resistivity is assumed to be constant on each patch. Equation (22) is plugged into the 


equation (55), and after integration, 





E,(P,, ) = —jkno 


~jkr N x y 
e {oe 
4nr 


a Ta 54 5 (gv) (57) 


n=) 


where the last factor in the equation is 





, Vv ° kAx . kA 
F (2 Vm | = oF ntim+ nvm) AxAry sine ste sine{ Am (58) 


ZZ 


Equation (57) can be solved for the unknown current coefficients given a set of m = 


1,2,.....N scattered field values. In matrix form, equation (57) is 


Fp F 


E; = E;=7-P (59) 


"by 


These two equations can be combined and put into the following form : 


EI /s Bf 7 (60) 


E; 
In equation (60), the only unknowns are the I* and I’ coefficients, and they can be found 
by simple matrix algebra. After these coefficients are found, they are plugged into the 


resistive sheet integral equation (equation (21)) to find the resistivity coefficients in 


1 _ g | he Ts J cr 
Vv’ 7 g I’ z ZX Zp (61) 


There are 4N rows in equation (61). In an isotropic resistivity case, there are 3 N 


equation (56). 


unknowns because the terms # “and # ” in equation (61) are same. Rearranging equation 


(61) 


gy LC’ V* Tx 7. i: 
, = = (62) 
Ce Tele Aye ZenieL” | I 
The matrix operations on the right hand side yield a column vector. The elements of the 


column vector are 77/7, and 77 J* and since I;, and I;, are known, the resistivity can 


mom 


be found by element-wise division. 
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IV. COMPUTER IMPLEMENTATION AND RESULTS 


A. INTRODUCTION 


In this chapter, the computer coding of the RCS synthesis formulas is described and 
numerical results are presented for a flat resistive sheet. The results of three investigations 


are included : 


1. Verification of the MM solution for rectangular patches and pulse basis 
functions. 


2. Verification of the synthesis equations derived in Chapter III and a parametric 
study of convergence. 


3. Sampling requirements. 


B. VALIDATION OF THE RECTANGULAR PATCH CODE 


The rectangular patch code based on pulse expansion functions developed in 
Chapter III was programmed in MATLAB. Computed results were compared to a MM 
code named PATCH [9]. Typical results are shown in Figures 6 through 10. The two codes 
give essentially the same RCS as expected. The agreement is within 2 dB down to the -40 
dB level when the two codes are converged. 

Validation of the rectangular patch code is important because many of the same 
subroutines and functions are used in the synthesis program. Furthermore, the output of this 
program is used as the input for initializing the synthesis process. Figures 11 and 12 show 
two dimensional plots of the RCS. The corresponding electric field values are used to fill 
the vector V in Equation (40). 


CG: VERIFICATION OF THE SYNTHESIS EQUATIONS 


The synthesis equations and computer codes are validated by the following process : 


1. Compute the scattered field for a plate with known surface resistivity. 
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2. Use the scattered field from step 1 in the synthesis procedure. 


3. Compare the synthesized resistivity in step 2 to the original resistivity used in 
step 1. In principle the difference should be zero. However, because of 
numerical and sampling errors, the residual will not be exactly zero, but some 
small (and hopefully negligible) value. 


The difference between the original and synthesized resistivity AR, is defined by 
A Ry ee 
where 
R,. = Onginal resistivity used to compute the electric field values. 


R,, = Synthesized resistivity. 


For R,,=377, AR, <107'* as shown in Figure 13. For R,=10, AR, is as large as 0.06 
as illustrated in Figure 14. Similar results are shown in Figures 15 through 19. The lower 
plate resistivity supports surface waves which increase the size of the off-diagonal elements 
in Z relative to the size of the diagonal elements. This leads to increased roundoff error. The 
number of field (RCS) samples NV, and AN, are equal, and they include the points outside 
of the visible region, which is the second approach as described in Chapter III-B. The 


number of patches are NV, = N, and N, = N,,. All of the visible region field points were 


in the backscatter (rear) hemisphere. Including the points outside the visible region affects 
the accuracy of the computation and causes some poorly conditioned matrices. However, 
the error resulting from this approach is not significant. In all cases the residual resistivity 


was less than 0.6 percent. 
D. SAMPLING REQUIREMENTS 


Computer memory and execution time are very important limitations in the 
scattering and synthesis computations. Reducing computer run time depends on reducing 
the size of the matrices (especially impedance matrix), and consequently minimizing the 
number of patches. Thus a guideline for the minkmum number of subdomains which gives 


computationally accurate results is necessary to minimize the computer run time. 


26 


For the synthesis calculation the continuous scattered field is sampled just as a 
continuous time function 1s sampled in a digital communication system. Therefore the same 
sampling requirements apply. If the frequency content of the scattered field is plotted, the 
minimum sampling rate must be at least twice the highest frequency [3]. To obtain the 
frequency content a two-dimensional discrete Fourier transform (DFT) is performed on the 
scattered electric field values, which are computed in sufficiently large number. The 
computation is performed in DCS because most RCS patterns are highly periodic in DCS 
and thus the number of sampling points is smaller. If N,, and Ny are the number of 
sampling points in the uw and v coordinates, then the sampling intervals are 2/N,, and 2/Ny. 
This represents 7 radians in visible space, and hence the scaling of the uw and v axes upon 
returning from the FFT is 2/N,, Auand 2/N,, Av, respectively. 

In Figures 20 through 23, frequency information for 1A x1A and 3A x 3A plates is 
plotted. The FFT is performed on the center row of the two-dimensional scattered electric 
field matrix to illustrate the frequency content in one coordinate. The magnitude of this 
transform is plotted with respect to the number of cycles in DCS in Figures 20 and 22. The 
result of a two-dimensional FFT performed on the whole E-field matrix is plotted with 
respect to the number of cycles in DCS in Figures 21 and 23. The maximum frequency for 
1X. x 1A resistive plate is 1 for the bistatic case. Since the sampling frequency must be twice 
the maximum frequency component of the signal, the sampling frequency is 2, which 
corresponds to 0.54. x 0.54 patch size. Applying this procedure to a 3A x 3A plate yielded 
the same maximum patch size. The required sampling frequency approaches zero as the 
plate size gets smaller, as expected. Thus in the case of bistatic RCS synthesis a half- 
wavelength subdomain size corresponds to the Nyquist sampling rate. Note that this 
guideline applies only to cases where the RCS pattern is highly periodic (1.e., negligible 
multiple reflections, surface waves, diffraction, etc.). When these other scattering 
mechanisms are present the scattering pattern becomes more complicated and thus contains 


more high-frequency components. Also, it is important to note that satisfying the field 


Zi 


sampling rate does not insure that the series expansion for the current on the plate is 


converged. 


Solid line - Patch code 





-30 -60 -40 = -20 0 20 40 60 80 
@, degrees (4 - 0°) 


Figure 6. Comparison of the bistatic RCS computed with PATCH vs that with 
rectangular patches (Rs = 377, Ny= Ny = 8, Ly= Ly= lA and G= gj= 0°). 
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Solid line - Patch code 





-30 -60 -40— -20 0 20 40 60 80 
@, degrees (4 = 0°) 


Figure 7. Comparison of the bistatic RCS computed with PATCH vs that with rectangular 
patches (Rs= 377, N,= N,= 14, L,= L= lA and @= ¢= 0°). 
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Solid line - Patch code 





-30 -60 -40-— -20 0 20 40 60 80 
@, degrees (4 = 0°) 


Figure 8. Comparison of the bistatic RCS computed with PATCH vs that with rectangular 
patches (R,= 10, N,= N= 20, L,=L,= lA and G= g= 0°). 
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Solid line - Patch code 


fs in dB 





-0 -60 -40 = -20 0 20 40 60 80 
@, degrees (¢ = 0°) 


Figure 9. Comparison of the bistatic RCS computed with PATCH vs that with rectangular 
patches (R,= 377, N= N,= 20, L,= Ly= lA and @= 30°, d= 0°). 
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Solid line - Patch code 





-30 -60 -40-— -20 0 20 840 60 SO 
@, degrees (4 = 0°) 


Figure 10. Comparison of the bistatic RCS computed with PATCH vs that with rectangular 
patches (K,= 377, N= N,= 20, L,=1,=2A,andiei= p= 0%). 
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Figure 11. Two dimensional plot of the bistatic RCS of a rectangular resistive sheet (R, 
Si ie = 23, (N= N,, NN) L.=L,= 1A and & = 


Q7- 0"). 
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Figure 12. Two dimensional plot of the bistatic RCS of a rectangular resistive sheet (R= 
377, Ly = Ly= lAand G= 30°, gj= 0°). 
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Figure 13. Residual resistivity for Rsg= 377, Ny=Ny=4 (Nyx=Ny, Ny=Ny), Ly = Ly= 1A and 
normal incidence. 
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Figure 14. Residual resistivity for Rsg= 10, Ny= Ny= 9 (N= Ny, Ny= Ny), Ly = Ly= 1a 
and normal incidence. 
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Figure 17. Residual resistivity for Rsg= 377, Ny= Ny= 16 (Ny= Ny, Ny= Ny), Ly = Ly= 2A 
and normal incidence. 
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Figure 18. Residual resistivity for Rsg= 377, Ny= Ny= 23 (Nyx= Ny, Ny= Ny), Ly = Ly= 1A 
ormal incidence. 
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Figure 20. Frequency analysis of one row (or uw coordinate) of the scattered E-field. Maximum 


1A, 256 field points and normal incidence. 
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Figure 22. Frequency analysis of one row (or one column) of the scattered E-field. Maximum 
frequency is 3 for R,.= 377, L,.= L,= 3A, 900 field points and normal incidence. 
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Figure 23. Frequency analysis in two-dimension for the scattered E-field. Maximum 


3A, 900 field points and normal 


frequency is 3 for Rsg= 377, Ly 


incidence. 
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V. CONCLUSIONS 


A rigorous bistatic RCS synthesis method for flat rectangular sheets has been 
presented. For convenience, the research was limited to resistive surfaces and parallel 
polarization. The scattering and synthesis equations are based on the method of moments 
(MM) technique. A rigorous method requires that both the current and resistivity be 
expanded into series, and the unknown coefficients determined simultaneously. 

The scattered field for a resistive sheet was first calculated for a known resistivity. 
The bistatic RCS of the resistive sheet was used to validate the rectangular patch solution 
by comparing it with the output of PATCH for different cases. This direct problem solution 
was used to generate and validate the synthesis procedure. The unknown current and 
resistivity coefficients were found by solving two coupled integral equations (equations (15) 
and (16) ). The reconstructed and onginal resistivities were compared to verify the 
synthesis equations. A parametric study of the solution convergence was performed. In all 
cases, the reconstructed resistivity obtained from the synthesis equations was in agreement 
with the original resistivity. 

The frequency content of the scattered electric field was also investigated in order to 
determine the maximum allowable patch size. The number of far-field points must be 
equal to the number of patches on the plate to obtain a unique solution. If the field pattern is 
highly periodic in DCS, as in the case of the physical optics solution, the Nyquist sampling 
theorem yields approximately 0.54 patch size on the plate for bistatic RCS. However, a 
MM solution with 0.54 patches is usually not converged; the conventional guideline is 
0.14 subdomains. However, in many cases of practical interest where the patterns to be 
synthesized are primarily determined by specular scattering (i.e., negligible surface waves 
and multiple reflection and diffraction) the minimum patch size can possibly be larger than 
the usual 0.12. 

The rigorous solution presented in this thesis is limited to flat rectangular resistive 


sheets. This allowed the use of rectangular subdomains with pulse basis functions. Pulse 
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basis functions are convenient for computer programming and result in matrix synthesis 
equations that are easy to solve. However, the more important scattering problems involve 
doubly curved surfaces. For these targets triangular subdomains must be used to accurately 
represent the surface. The overlapping basis functions of the type used in PATCH will 
result in much more complicated (and possibly nonlinear) sets of equations. It is 


recommended that this be the direction of future research in this area. 
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APPENDIX COMPUTER CODES 


In this appendix computer codes are presented. The code plate.m calculates the 
RCS in @-direction for the bistatic case and TM, polarization. The code platedcs.m 
calculates the RCS and plots it in direction cosine space. The data from this code is saved 
to be used in the synthesis program. The code synt.m uses the stored data to reconstruct the 
surface resistivity. The codes dcsin.m, center.m, green.m, greenl.m, green2.m, green3.m, 
and green4.m are the functions which are used in the synthesis and scattering programs. 
The code sample.m performs the FFT analysis in the u and v directions of the scattered 


field. Detailed explanations are presented in the comments in the codes. 
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WALWPVA/A///PAVPAV/A’/VVAD/VVVVVVVVAVWAVVVMU%%% 


% plate.m % 
% Resistive Square Plate Bistatic RCS (in teta direction) % 
% with Method of Moments. % 
% % 
% Ugurcan Samli, July 1996 % 
% % 


MPAPPP///// PLP P/V VV % VV A/VVWAVAVAVAV%V% 


clear all; 
format compact; 


PWV AHVVWVVVWVV%V%% Variables used in the program%%%%%%%%%%%%% 


’% modex and modey are the number of patches in x- and y-directions 
% respectively. Eotet and Eophi are for polarization. Rs is the 
% constant surface resistivity. 


Lx=1; % input(‘Length in x-direction= '); 
Ly=1; % input(‘Length in y-direction= '); 
modex=20; modey=20; 
deltax=Lx/modex; deltay=Ly/modey; 
mode=modex*modey; eta=377; lambda=1; 

=2* pi; t=pi/180; Rreal=377; Rimg=0; 
Eotet=1; Eophi=0; Emag=sqrt(Eotet*2+Eophi‘2); 
incr=2; sta=-90; stp=90; 
lim=floor((stp-stat+1e-5)/incr)+1; 
int=2; phiang=0; phi=phiang*t; 


Rsl=Rreal+j*Rimg; Rs2=Rs1*eye(mode); Rs=diag(Rs2)'; 
PAWVVMV%%%Y% Computing patch center points %%%%%% 
in=0; 

for s=1:modey 

for ss=1:modex 


in=int1; 


xn(in)=(-Lx/2+deltax/2+(ss-1)*deltax); 
yn(in)=(-Ly/2+deltay/2+(s-1)*deltay); 


end; 
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end; 

xm=xn; ym=yn; 

funl=j*k*eta; fun2=j/k* eta; 

PPPVYV/%V%Y Calculation of impedance matrix %%%%%%%%%%% 


% Impedance matrix calculation is same with the impedance matrix 
% calculation in platedcs.m, green functions are used in calculations. 


for n=1:mode 
for m=1:mode 


Gxx=green(xm(m),xn(n),ym(m),yn(n),deltax,deltay,int); 
Gyy=green(ym(m),yn(n),xm(m),xn(n),deltay,deltax,int); 


xn1=xn(n)+deltax/2; 
xn2=xn(n)-deltax/2; 
xm1=xm(m)+deltax/2; 
xm2=xm(m)-deltax/2; 
ynl=yn(n)+deltay/2; 
yn2=yn(n)-deltay/2; 
ym1=ym(m)+deltay/2; 
ym2=ym(m)-deltay/2; 


Qxx=(green1(xm2,xn2,ym(m),yn(n),deltax,deltay,int)... 
-greenl (xm2,xn1,ym(m),yn(n),deltax,deltay,int)... 
-greenl(xm1,xn2,ym(m),yn(n),deltax,deltay,int)... 

+green1(xm1,xnl,ym(m),yn(n),deltax,deltay,int)); 

Oxy=(green2(xm2,xn(n),ym(m),yn2,deltax,deltay, int)... 
-green2(xm2,xn(n),ym(m),yn1 ,deltax,deltay,int)... 
-green2(xm1,xn(n),ym(m),yn2,deltax,deltay,int)... 
+green2(xm1,xn(n),ym(m),yn1,deltax,deltay,int)); 

Qyx=(green3(xm(m),xn2,ym2, yn(n),deltax,deltay,int)... 
-green3(xm(m),xn2,ym1,yn(n),deltax,deltay,int)... 
-green3(xm(m),xn1,ym2,yn(n),deltax,deltay,int)... 

+sreen3(xm(m),xnl,ym1,yn(n),deltax,deltay,int)); 

Qyy=(green4(xm(m),xn(n),ym2,yn2,deltax,deltay,int)... 

-green4(xm(m),xn(n),ym2,yn1 ,deltax,deltay,int)... 
-green4(xm(m),xn(n),ym1,yn2,deltax,deltay,int)... 
+green4(xm(m),xn(n),ym1,yn1,deltax,deltay,int)); 


5] 


sum=fun 1] *Gxx-fun2*Qxx; 
if m==n, sum=sum+Rs(m)*(deltax/deltay); end 
Z(m,n)=sum; 


sum=-fun2*Qxy; 
Z(m,n+mode)=sum; 


sum=-fun2* Qyx; 
Z(m+mode,n)=sum; 


sum=fun1*Gyy-fun2*Qyy; 
if m+mode—=n+mode, sum=sum+Rs(m)*(deltay/deltax); end 
Z(m+mode,nt+mode)=sum; 


end; 
end; 


Zl=pinv(Z); 
disp(‘Impedence Matrix Computed and Inverted’); 


WWV/PV//VWVV/%Yo Incidence angle %%%%%%%%%%% 


phii=0; teti=0; 
sinphii=sin(phi*t); sinteti=sin(teti*t); 
cosphii=cos(phii*t); costeti=cos(teti*t); 


vvi=sinphii* sintet; 
uui=sinteti* cosphi1; 


Ex=Eotet* costeti*cosphii-Eophi* sinphii; 
Ey=Eotet*costeti* sinphiit+Eophi*cosphil; 


for n=1:mode; 
argex=k*deltax*uui/2; 
argy=k*deltay*vvi/2; 
sincx=1; sincy=1; 
if abs(argx)>1e-5, sincx=sin(argx)/argx; end 
if abs(argy)>le-5, sincy=sin(argy)/argy; end 
expxy=exp(j*k*(xn(n)*uuityn(n)*vvi)); 
vx(n)=Ex*deltax* expxy *sincx* sincy; 
vy(n)=Ey *deltay*expxy*sincx* sincy; 


a2 


end; 
% Excitation vector V, and current coefficients I. 


V(1:mode)=vx(:); 
V(mode+1:2*mode)=vy(:); 
I=ZI* V..'; 


%)%4%%H>%H%M%%™% calculation of the scattered fields %%%%%%%%% 


for aa=1:lim; 
angle(aa)=(aa-1)*incr+sta; 
teta=angle(aa)*t; 
costet=cos(teta); 
sintet=sin(teta); 
cosphi=cos(phi); 
sinphi=sin(ph1); 
vv=sinphi* sintet; 
uu=sintet* cosphi; 


for n=1:mode 
argx=k* deltax*uu/2; 
argy=k* deltay*vv/2; 
sincx=1; sincy=1; 
if abs(argx)>le-5, sincx=sin(argx)/argx; end 
if abs(argy)>1le-5, sincy=sin(argy)/argy; end 
expxy=exp(j*k*(xn(n)*uutyn(n)*vv)); 
vx 1(n)=deltax* expxy*sincx* sincy; 
vy 1(n)=deltay*expxy*sincx*sincy; 
end; 


Exs=vx1*I(1:mode); 
Eys=vy 1*I(mode+1:2*mode); 


%%%% TM and TE polarization cases are together.%%%%% 
etet=-)*k*eta/2/k*(Exs*costet* cosphitEys*costet* sinphi); 
ephi=-j*k*eta/2/k*(-Exs*sinphit+Eys*cosphi); 
etetrcs(aa)=etet; 


ephircs(aa)=ephi; 


end; 
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HV%%%%% Calculation of the RCS, and plots %%%%%%% 


RCStet=10* log10(abs(etetrcs).“2*2*k+1e-10); 
RCSphi=10*log10(abs(ephircs).*2*2*k+1e-10); 


plot(angle,RCStet), grid; 

max 1=max(max(RCStet)); 
db=(floor(max(max1)/5)+1)*5; 
axis([-90,90,db-50,db]); 

xlabel('‘Angle in degrees’); 

ylabel((RCS in dBsm ‘); 

title(‘Bistatic RCS of a Resistive Plate With MM’); 


Ysave f1 Z RCStet angle Rreal 
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LAW V/V/V V/V/V AYALA AY/AYS AYA Ve 


% platedcs.m 

% This program calculates Resistive Square Plate Bistatic RCS with 
% Moment Method, and it plots this RCS in direction cosine space. 
% 

% Ugurcan Samll. 

% 


“% 
% 
% 
% 
% 
% 


LUWWLLW% VAAL ALAA AY ALAA ASA 


clear all; 
format compact; 


WMWWVV%VV%%%V%Y Variables used in the program %%%%%%%%%%% 


IEx—1; % input(‘Length in x-direction='); 
Ly=1; % input(‘Length in y-direction= '); 
T=1; % If the number of field points equal to number of 


% patches, T=1. 
% modex and modey are the number of segments in x and y direction. 
% Eotet and Eophi determines the polarization (Eotet=1 and Eophi=0 
% is the Teta polarization.) 
eta=377; lambda=1; k=2* pi; t=pi/180; int=2; 
Rreal=377; Rimg=0; 
Eotet=1; Eophi=0; 


% Finding center of the patches in Direction Cosine Space. 
% dcsin.m and center.m functions are called. 

modex=4; modey=3; 

deltax=Lx/modex; deltay=Ly/modey; 

mode=modex* modey; 

load uv1l % u and v points are computed seperately. 


% Constant surface resistivity Rs matrix. 
Rs1=Rreal+j*Rimg; Rs2=Rs1*eye(mode); Rs=diag(Rs2)'; 


WWMHVVVVH%%% Computing Patch center points. %%%%%%%%%%%%% 


in=0; 

for s=1:modey 

for ss=1:modex 
in=in+1; 
xn(in)=(-Lx/2+deltax/2+(ss-1)*deltax); 


>> 


yn(in)=(-Ly/2+deltay/2+(s-1)*deltay); 
end; 
end; 


xm=xn; ym=yn; 
funl=j*k* eta; fun2=)/k*eta; 
MMV%%%% Filling Impedence Matrix %%%%%%%%%%M%% 


% Green functions are called in the impedance matrix calculation. 
for n=1:mode 
for m=1:mode 


Gxx=green(xm(m) ,xn(n) sym(m) syn(n) ,deltax,deltay, int) : 
Gyy=green(ym(m),yn(n),xm(m),xn(n),deltay,deltax, int); 


xn1=xn(n)+deltax/2; 
xn2=xn(n)-deltax/2; 
xm 1=xm(m)+deltax/2; 
xm2=xm(m)-deltax/2; 
ynl=yn(n)+deltay/2; 
yn2=yn(n)-deltay/2; 
ym1l=ym(m)+deltay/2; 
ym2=ym(m)-deltay/2; 


Qxx=(green1(xm2,xn2,ym(m),yn(n),deltax,deltay,int)... 
-greenl(xm2,xn1,ym(m),yn(n),deltax,deltay,int)... 
-greenl(xm1,xn2,ym(m),yn(n),deltax,deltay,int)... 
+greenl(xm1,xnl,ym(m),yn(n),deltax,deltay,int)); 

Qxy=(green2(xm2,xn(n),ym(m),yn2,deltax,deltay,int)... 
-green2(xm2,xn(n),ym(m),yn1 ,deltax,deltay,int)... 
-green2(xm1,xn(n),ym(m),yn2,deltax,deltay,int)... 
+green2(xm1,xn(n),ym(m),yn1 ,deltax,deltay,int)); 

Qyx=(green3(xm(m),xn2,ym2, yn(n),deltax,deltay,int)... 
-green3(xm(m),xn2,ym1,yn(n),deltax,deltay,int)... 
-green3(xm(m),xnl,ym2,yn(n),deltax,deltay,int)... 
+green3(xm(m),xnl,ym1,yn(n),deltax,deltay,int)); 

Qyy=(green4(xm(m),xn(n),ym2,yn2,deltax,deltay,int)... 
-green4(xm(m),xn(n),ym2,yn1,deltax,deltay,int)... 
-green4(xm(m),xn(n),ym1,yn2,deltax,deltay,int)... 
+green4(xm(m),xn(n),ym1,yn1,deltax,deltay,int)); 
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end; 


sum=fun] * Gxx-fun2*Qxx; 
if m==n, sum=sum+Rs(m)*(deltax/deltay); end 
Z(m,n)=sum; 


sum=-fun2* Qxy; 
Z(m,n+mode)=sum; 


sum=-fun2*Qyx; 
Z(m+mode,n)=sum; 


sum=fun1*Gyy-fun2*Qyy; 
if m+mode==n+mode, sum=sum+Rs(m)*(deltay/deltax); end 


Z(m+mode,n+mode)=sum; 


end; 


Zi=inv(Z); 
disp(Impedence Matrix Computed and Inverted’); 


MMV%%%% Incident Angles %%%%w%%%% 


phii=0; teti=0; 
sinphii=sin(phii*t); sinteti=sin(teti*t); 
cosphi=cos(phu*t); costeti=cos(teti*t); 


vvi=sinphii* sinteti; 
uui=sinteti* cosphil; 


Ex=Eotet* costeti*cosphiu-Eophi*sinphii; 
Ey=Eotet* costeti* sinphiit+Eophi*cosphu; 


for n=1:mode; 
argx=k* deltax*uur/2; 
argy=k* deltay*vvi/2; 
sincx=1; sincy=1; 
if abs(argx)>1e-5, sincx=sin(argx)/argx; end; 
if abs(argy)>1le-5, sincy=sin(argy)/argy; end; 
expxy=exp(J*k*(xn(n)*uultyn(n)*wv1)); 
vx(n)=Ex* deltax* expxy*sincx* sincy; 
vy(n)=Ey*deltay* expxy*sincx* sincy; 


Sy 


end; 


% Excitation vector and current coefficients, V and I respectively. 
V(1:mode)=vx(:); 

V(modet1 :2*mode)=vy(:); 

[=Z1*V.'; 


WV %%% calculation of the field points%%%%%%%% 


for m=1:mode 

unit 1 (m)=(sqrt(u(m)2+v(m)2)); 
phi(m)=(atan2(u(m)+1e-7,v(m))); 
teta(m)=asin(unit1(m)); 
costet=cos((teta(m))); 
sintet=sin((teta(m))); 
cosphi=cos(phi(m)); 
sinphi=sin(phi(m)); 


for n=1:mode 

argx=k* deltax*u(m)/2; 

argy=k* deltay* v(m)/2; 

sincx=1; sincy=1; 

if abs(argx)>1e-5, sincx=sin(argx)/argx; end 
if abs(argy)>1le-5, sincy=sin(argy)/argy; end 
expxy=exp(j*k*(xn(n)*u(m)+yn(n)*v(m))); 
vx 1(n)=deltax*expxy*sincx* sincy; 

vy 1(n)=deltay*expxy*sincx* sincy; 
end; 


% x- and y-components of the scattered field. 
Exs=vx1*I(1:mode); Exs1(m)=Exs; 


Eys=vy 1 *I(mode+1:2*mode); Eysl(m)=Eys; 


etet=-]*k* eta/2/k*(Exs*costet*cosphitEys*costet*sinphi); etet](m)=etet; 
ephi=-)*k* eta/2/k*(-Exs*sinphit+Eys*cosphi); ephil(m)=ephi; 


etetrcs(m)=(-]*k*eta/2/k)*Exs; 
ephircs(m)=(-]*k*eta/2/k)*(Eys); 


end; 


~MAVHVWVV% Calculation of RCS %%%%%%%%% 
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for i=l:mode 
RCStet(i)=10*log10(abs(etetres(i)).~2*2*k+1e-10); 
RCSphi(i)=10* log 10(abs(ephircs(1)).“2*2*k+1e-10); 
end; 

RCS=reshape(RCStet,modex,modey); 
u=reshape(u,modex,modey); v=reshape(v,modex,modey); 


mesh(u,v,RCS); 


% All the data are saved for use in synthesis program. 
save b24 RCS u v mode modex modey Exs] Eys1 V Z deltax deltay xn xm yn ym Rs I 
etetrcs ephircs etet] ephil etet ephi phii teti RCSphi 
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WLP V0%0%0 700% 0/00 V0 VOY 0VOV0VOVOVOVO YOO VVV%VVVVVV%V%V%0% 


% synt.m 

% 

% This program calculates the synthesis results for bistatic case. 
% The result will be the resistivity which is found 

% from the synthesis equations. 

% 

% Ugurcan Samli, July 1996 

% 


Yo 
% 
Yo 
% 
% 
% 
% 
% 


WWVVVV/VOY% V0 V0 P0000 VV VOVOYO OVO VO VOY VV VOVVVVV%%V0%Y% 


clear all; 
format long; 


% Reading in the data file from the scatterring program (platedcs.m) 


load b24.mat; 
disp('Data file is loaded.'); 


eta=377; lambda=1; t=p1/180; k=2* pi; 
uxx=linspace(-1,1,mode); 

vyy=linspace(-1,1,mode); 

% Reading the E-field components from the data file. 


efld(1:mode)=Exs1(:); efld(mode+1:2*mode)=Eys1(:); 
rs=[Rs,Rs]; rs=diag(rs); 


PWVWVVV%V%%%% Filling the matrices %%%%%%%%%%%%% 
% Oniginal resistivity is subtracted from the impedance matix. 
Z=Z-IS; 
for m=1:mode 

for n=1:mode 

argy=k* deltay*v(m)/2; 


argx=k*deltax*u(m)/2; 
sincx=1; sincy=1; 
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if abs(argx)>1e-5, sincx=sin(argx)/argx; end; 
if abs(argy)>1e-5, sincy=sin(argy)/argy; end; 
expxy=exp(j*k*(xn(n)*u(m)+yn(n)*v(m))); 


fnm 1(m,n)=expxy*sincx*sincy*deltax; 
fnm2(m,n)=expxy*sincx* sincy *deltay; 


end; 
end; 


% Filling the scattered matrix (First equation), and current 
% coefficients, I. 


FF 1=pinv(fnm 1); FF2=pinv(fnm_2); 
FFI=[FF 1,zeros(mode,mode);zeros(mode,mode),FF2]; 


I=(efld*FFI.'); 
MWHVMNM%V%%%%% Matrix Operations %%%%%%%%%%% 


% Second equation (Equation 63) operations. Excitation vector V is from the 
% data file. 


Zixy=Z*1.'; 
AA=V.'-ZIxy; 
ZL=(AA,/L'); 
ZL1=real(ZL); 


ZLx=(deltay/deltax).*ZL1(1:mode); 
ZLy=(deltax/deltay).*ZL 1(mode+1 :2* mode); 


% Reconstructed resistivity rx. 
rx=reshape(ZLx,modex,modey); 
rorg=reshape(Rs,modex,modey); 


IX=1X-rorg; 


uxx=reshape(uxx,modey,modex); vyy=reshape(vyy,modex,modey); 
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WMV AVY Plots *VAVAVAVMVAVAVAVAV%%% 
% Plot of residue (rx) resistivity. 

mesh(uxx',vyy,rx); grid; view(30,30); 

“save s12 uxx vyy rx rorg 


title’Comparision of the Resistivities; in x-direction’); 
xlabel(‘u’), ylabel('v'), zlabel('Resistivity Differences’); 
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function [nu,nv|=dcsin(mode) 


~AKUWAWWVWWWVV/VV/YV%/ALVV/VWA/VWVWV%/AV%//%N%Y/0% 


% This function calculates the number of points inside the 


% unit circle, and shows these points. 
% 

% function [nu,nv]=dcsin(mode) 

% 

% Ugurcan Samli 


% 
% 
% 
% 
% 
% 


UWA V/V/V V/V/V /%%0%0% 


% Area ratio of the unit circle and the total area. 


tot_pts=1.28*mode; 
totsqr=sart(tot_pts); 


if abs(round(totsqr)-totsqr)<0.2 
nu=round(totsqr); 
nv=round(totsqr); 

else 
nu=ceil(totsqr); 
nv=floor(totsqr); 

end; 


[u,v]=center(nu,nv); 


n uv=nu* nv; 


sintet=sqrt(u.2+v.%2); 
for s=l:n_uv 


if sintet(s)>1 
u(s)=0; 
v(s)=0; 
else 
u(s)=u(s); 
v(s)=v(s); 
end; 


end; 


aa=find(u=0); bb=find(v—=0); 


63 


u(aa)=[]; v(bb)=[]; 
[nu nv length(u) length(v)} 


plot(u,v,'+'),xlabel(‘u’),ylabel(‘v'); axis([-1,1,-1,1]); 
title(‘Points inside the unit circle’); 
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function [x,y |=center(modex,modey) 


WVWVVVV/ VV V0VOVOVOVO VO VO VOYVO%VO%%YOV0VOVOVOV0 00 VO VOY%OVO%OY0 


% center.m % 
% This Function finds the center points of the patches % 
% in the Direction Cosine Space. % 
% % 
% [x,y ]=center(modex,modey) 

% 

% % 
% Ugurcan Samli, May, 1996. % 


WAMWAWVAWVVWVVVVWVAVVAVVVVVVVAVVVAVAVVVMUVVVV%%0% 
Lx=2; Ly=2; 

deltax=Lx/modex; deltay=Ly/modey; 

a=0; 


for s=1:modey 
for ss=1:modex 


a=a+1; 
x(a)=(-Lx/2+deltax/4+(ss-1)*deltax); 
y(a)=(-Ly/2+deltay/4+(s-1)*deltay); 


end; 
end; 
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Y% 


green.m 


% Ugurcan Samli, Resistive Plate MM solution 
% This function computes the integral of green's function. 


function gg=green(a,b,c,d,g,h,int); 


w—0.225*h: 


if int==0 


end 


r=sqrt(((b-a)*2+(c-d)“2)+w’2); 
ggl=exp(-2*)*pi*r)/r; 
ge=gol*9/4/pi*g; 


if int—=2 


end 


xg=[0.5773503 -0.5773503]; 
ag(1:2)=[1,1]; 


if int~=0 


int2=1nt/2; 
for s=] :int2 
sl=int-s+1; 
ag(sl)=ag(s); 
xg(sl)=xg(s); 
xg(s)=-xg(s); 


end; 


end; 


xhigh=a+¢/2; 
xlow=a-¢/2; 
bb=(xhigh-xlow)/2; 
aa=(xhightxlow)/2; 
sum=0; 


for 1=]:1nt 
x=bb* xg(i)+aa; 
r=sqrt((b-x)2+(c-d)\2+w%2); 
sum=sum+ag(1)*exp(-j*2* pi*r)/r; 
end; 
gg=sum*bb/4/pi*g; 
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% greenl.m 
% Ugurcan Samli, Resistive Plate MM solution 
% This function computes the integral of green's function. 


function gg=greenl(a,b,c,d,g,h,int); 
w=0.225*¢; 


iint——0 
r=sqrt(((a-b)*2+(c-d)*2)+w’2); 
ggl=exp(-2*)*pi*r)/r; 
gg=ge1*h/4/pi/h; 

end; 


if int—=2 
xeg=[0.5773503 -0.5773503]; 
ag(1:2)=[1,1]; 

end; 


if int~=0 
yhigh=ct+h/2; 
ylow=c-h/2; 
bb=(yhigh-ylow)/2; 
=(yhight+ylow)/2; 
sum=0; 
for 1=1:int 
y=bb*xg(i)+aa; 
r=sqrt((a-b)*2+(y-d)\2+w"2); 
sum=sum-+ag(1)*exp(-j*2*pi*r)/r; 
end; 
ge=sum*bb/4/pi/h; 
end; 
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% 


green2.m 


% Ugurcan Samli, Resistive Plate MM solution 
% This function computes the integral of green's function. 


function Bom een2(a,b,c,d, g,h,int); 


w=0.225* 9; 


if int—0 


end; 


r=saqrt(((a-b)*2+(c-d)*2)+w’2); 
ggl=exp(-2*)*pi*r)/r; 
ge=g¢gl*h/4/pi/h; 


ie 


end; 


xg=[0.5773503 -0.5773503]; 
ag(1:2)=[1,1]; 


if int~=0 


end; 


int2=int/2; 
for s=1:int2 
sl=int-s+1; 
ag(s1)=ag(s); 
xg(s1)=xg(s); 
xg(s)=-xg(s); 
end; 
yhigh=c+h/2; 
ylow=c-h/2; 
bb=(yhigh-ylow)/2; 
aa=(yhight+ylow)/2; 
sum=0; 
for 1=1] :int 
y=bb*xg(i)+aa; 
r=sqrt((a-b)*2+(y-d)*2+w%2); 
sum=sumtag(1)*exp(-j*2*pi* rr; 
end; 
ge=sum*bb/4/pi/h; 
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% 


green3.m 


% Ugurcan Samli Resistive Plate MM solution 
% This function computes the integral of green's function. 


function ge=green3(a,b,c,d,g,h,int); 


w=(0.225*h; 


if int==0 


end; 


r=sqrt(((a-b)*2+(c-d)*2)+w”2); 
ggl=exp(-2*)*pi*r)/r; 
ge=-g2e1*g/4/pi/g; 


iint—2 


ena: 


xg=[0.5773503 -0.5773503]; 
ag(1:2)=[1,1]; 


if int~=0 


end; 


int2=int/2; 
for s=1 :int2 
sl=int-s+1; 
ag(s1)=ag(s); 
xg(s1)=xg(s); 
xg(s)=-xg(s); 
end; 
xhigh=at+¢/2; 
xlow=a-¢/2; 
bb=(xhigh-xlow)/2; 
=(xhightxlow)/2; 
sum=0; 
for 1=1:int 
x=bb*xg(i)t+aa; 
r=sqrt((x-b)\2+(c-d)\2+w’2); 
sum=sumtag(1)*exp(-j*2*pi*r)/r; 
end: 
ge=sum*bb/4/pi/g; 
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% 


green4.m 


% Ugurcan Samli Resistive Plate MM solution 
% This function computes the integral of green's function. 


function gg=green4(a,b,c,d,g,h, int); 


w=0.225*h; 


if int==0 


end; 


r=sqrt(((a-b)2+(c-d)*2)+w’2); 
go 1=exp(-2*)*pi*r)/r; 
ge=ge1*g/4/pi/g; 


if int==2 


end; 


xg=[0.5773503 -0.5773503]; 
ag(1:2)=[1,1]; 


if int~=0 


end; 


int2=int/2; 
for s=1:int2 
sl=int-s+1; 
ag(s1)=ag(s); 
xg(sl)=xg(s); 
xg(s)=-xg(s); 
end; 
xhigh=a+¢/2; 
xlow=a-g/2; 
bb=(xhigh-xlow)/2; 
aa=(xhight+xlow)/2; 
sum=0; 
for 1=1:1nt 
x=bb*xg(i)+aa; 
r=sqrt((x-b)“2+(c-d)*2+w%2); 
sum=sum+ag(1)*exp(-j*2* pi*r)/r; 
end; 
gg=sum*bb/4/pi/g; 
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PWVVVVWVVVV/0V%AV%07%/0%V%0%YO/0 VV VV %%%0%%0%0% 


% 

% sample.m 

% This program analyses the frequency content of the E-field 
%  inuand v directions. 

% 

%  Ugurcan Samli, August 96 

% 


% 
Yo 
% 
% 
% 
% 
% 


~LLY/AWWLVVV/VV%//%%%/ ALLA WLLVL/V%/7/A0%%% 


clear all; 
load b24.mat 


Efld=Exs1; 
Efld=reshape(Efld,modex1 ,modey 1); 


% fn : number of FFT points. k : Spatial Frequency index. 
% us : Sampling interval in DCS. 
% This is for bistatic case and normal incidence. 


fn=128; 
cl=0:fn-1; 


us=2/(modex1-1); 
k=c1*2/(fn*us); 


% fft in 1-dimension. FFT is performed only in 1 column (or row). 


k=k(1 :fn/2+1); 

el=Efld(6,:); 

e2=abs(fft(e1,fn)); 

figure(1); 

plot(k,e2(1:fn/2+1)); grid; 
xlabel(‘Number of samples in DCS’); 
ylabel(‘abs(fft)’); title(3 by 3 plate’); 


% fft in 2-dimension. 
ffabs(fft2(Efld,fn,fn)); 
[ku,kv]=meshgrid(k,k); 


figure(2); 
mesh(ku, kv, ff(1:fn/2+1,1:fn/2+1)), grid; view(30,30); 
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xlabel(’ Number of Samples in u-direction’); 
ylabel(‘Number of Samples in v-direction’), 
Zlabel(‘abs(fft)'); title(3 by 3 plate’); 


72 


LIST OF REFERENCES 


Faros, N. I., “Radar Cross Section Synthesis for Planar Resistive Surfaces,” 
Master’s Thesis, Naval Postgraduate School, 1994. 


Waddel, C. W., “Radar Cross Section Synthesis of Doubly Curved Surfaces,” 
Master’s Thesis, Naval Postgraduate School, 1995. 


David C. Jenn, Radar and Laser Cross Section Engineering, AIAA Education 
Series,Ohio, 1995. 


Balanis, A. C., Advanced Engineering Electromagnetics, John Wiley & Sons, 
Canada, 1989. 


“Survey of Numerical Electromagnetic Modeling Techniques,’ WWW Internet, 
University of Missouri, EMC Laboratory. 


Randy Bancroft, Understanding Electromagnetic Scattering Using the Moment 
Method, Artech House, Boston, 1996. 


D. C. Jenn, unpublished notes, Naval Postgraduate School, Monterey, CA. 


J. Moore and R. Pizer, Moment Methods in Electromagnetics, Research Studies 
Press, John Wiley & Sons Inc., June 1986. 


es 





INITIAL DISTRIBUTION LIST 


Defense Technical Information Center 
8725 John J. Kingman Rd., STE 0944 
Ft. Belvoir, VA 22060-6218 


Dudley Knox Library 

Naval Postgraduate School 

411 Dyer Rd. 

Monterey, California 93943-5101 


Professor David C. Jenn, Code EC/IN 

Department of Electrical and Computer Engineering 
Naval Postgraduate School 

Monterey, California 93943-5121 


Professor D. Curtis Schleher, Code EC/SC 
Department of Electrical and Computer Engineering 
Naval Postgraduate School 

Monterey, California 93943-5121 


Chairman, IW Academic Group 
Frederic H. Levien, Code EC/LV 
Naval Postgraduate School 
Monterey, California 93943 


Ltjg Ugurcan Samli 

Akmescit Sok. Ferah Apt. No: 1-11 
Bakirkéy-Istanbul 

Turkey 


Deniz Harp Okulu Komutanl1§1 Kiitiiphanesi 
Dz.H.O. Komutanlig1, Tuzla - Istanbul, 81704 
Turkey 


Yildiziar Suiistii Etit. Mrk. Kom.’l181 Elektronik Harp 
Destek Merkezi Merkezi Kom. ’lig1 

Golcuk - Kocaeli 

Turkey 


75 


No. Copies 
z 





DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY CA 93943-5101 





3 2768 00324493 07 


wei ii i 7 





